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Abstract 

The decay —>■ ’4}{2S)K~^Tr~ is observed using a data set corresponding to an 

integrated luminosity of 3.0 fb~^ collected by the LHCb experiment in pp collisions 
at centre-of-mass energies of 7 and 8 TeV. The branching fraction relative to the 
> 'ip{2S)K^TT~ decay mode is measured to be 

B(m Z = 5-38 ± 0.36 (stat) ± 0,22 (ayat) ±0.31 {/,//a) %, 

where fs/fd indicates the uncertainty due to the ratio of probabilities for a b 
quark to hadronise into a B^ or B^ meson. Using an amplitude analysis, the 
fraction of decays proceeding via an intermediate A* (892)*^ meson is measured to 
be 0.645 ± 0.049 (stat) ± 0.049 (syst) and its longitudinal polarisation fraction is 
0.524±0.056 (stat) ±0.029 (syst). The relative branching fraction for this component 
is determined to be 

B{B^ ^il;{2S)K*{Smf) ^ ^ ^ ^ 

H(R0 ^ V'(25')A*(892)0) ^ ^ ^^ \Js/jdj 

In addition, the mass splitting between the B^ and B^ mesons is measured as 

M(R°) - M(H°) = 87.45 ± 0.44 (stat) ± 0.09 (syst) MeV/c^. 
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1 Introduction 


The large data set collected by the LHCb experiment has allowed precision measurements of 
time-dependent CP violation in the —>■ J/'00 and 5° —)■ J/'0/o(98O) decay modes (HQ 

The results are interpreted assuming that these decays are dominated by colour-suppressed 
tree-level amplitudes (Fig. [^. Higher-order penguin amplitudes, which are difficult 
to calculate in QCD, also contribute (Fig. [^. Reference snggests that the size of 
contributions from these processes can be determined by studying decay modes snch as 
—)■ J/'^iF*(892)° where they dominate. The —)■ (892)° decay mode was first 

observed by the CDF collaboration and subseqnently stndied in detail by the LHCb 
collaboration [^. 



Figure 1: Tree (left) and penguin (right) topologies contributing to the —)• V’F decays where 

tP = J/P;^Pj{2S) and V = (/>,K*(892f. 


Recently, interest in 6-hadron decays to final states containing charmonia has been 
generated by the observation of the Z(4430)“ —)■ 'ijj{2S)7i~ state in the —)■ 'il>{2S)K~^7i~ 
decay chain by the Belle j6]-[^ and LHCb collaborations |^. As this state is charged 
and has a minimal qnark content of ccdu, it is interpreted as evidence for the existence 
of non-gg mesons 10 . Evidence for similar exotic strnctures in —)■ Xci,c2A'’''7r“ and 
R° —)■ J/ipK~^7r~decai'ys has been reported by the Belle collaboration 11,12. If these 


strnctnres correspond to real particles they shonld be visible in other decay modes. 

This letter reports the first observation of the decay B^ —>■ 'ip{2S)K~^7i~ and presents 
measurements of the inclusive branching fraction and the fraction of decays that proceed 
via an intermediate A*(892)° resonance, as determined from an amplitnde analysis of 
the final state. The amplitnde analysis also allows the determination of the longitndinal 
polarisation fraction of the A*(892)° meson. Additionally a measnrement of the mass 
difference between R° and B^ mesons is reported that improves the current knowledge of 
this observable. 


2 Detector and simulation 

The LHCb detector [I3p4| is a single-arm forward spectrometer covering the psendorapidity 
range 2 < r/ < 5, designed for the study of particles containing 6 or c quarks. The detector 

^Charge-conjugatation is implicit unless stated otherwise. 
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includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes 15 placed downstream of the magnet. The tracking system 
provides a measurement of momentum, p, of charged particles with a relative uncertainty 
that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The minimum distance 
of a track to a primary vertex, the impact parameter, is measured with a resolution of 
(15 -|- 29/pt) M-ni, where px is the component of the momentum transverse to the beam, 
in GeV/c. Large samples of —)■ and J/ijj decays, collected concurrently 

with the data set used here, were used to calibrate the momentum scale of the spectrometer 
to a precision of 0.03 ' 


i% 16 


Different types of charged hadrons are distinguished using information from two ring¬ 
imaging Cherenkov detectors 17 . Photons, electrons and hadrons are identihed by a 


calorimeter system consisting of scintillating-pad and preshower detectors, an electromag¬ 
netic calorimeter and a hadronic calorimeter. Muons are identihed by a system composed 
of alternating layers of iron and multiwire proportional chambers (^. The online event 
selection is performed by a trigger j^, which consists of a hardware stage, based on 
information from the calorimeter and muon systems, followed by a software stage, which 
applies a full event reconstruction. In this analysis candidates are hrst required to pass 
the hardware trigger, which selects muons and dimuon pairs based on the transverse 
momentum. At the subsequent software stage, events are triggered by a '0(2*5') —)■ 
candidate where the '0(25') is required to be consistent with coming from the decay of a 6 
hadron by either using impact parameter requirements on daughter tracks or detachment 
of the 0(25') candidate from the primary vertex. 

The analysis is performed using data corresponding to an integrated luminosity of 
1.0 fb“^ collected in pp collisions at a centre-of-mass energy of 7 TeV and 2.0 fb“^ collected 
at 8 TeV. In the simulation, pp collisions are generated using Pythia pO] with a specihc 


LHGb conhguration 21 . Decays of hadronic particles are described by EvtGen 22 


in which hnal state radiation is generated using Photos [^. The interaction of the 
generated particles with the detector and its response are implemented using the Geant4 
toolkit 


24 as described in Ref. 25 


3 Event selection 

The selection of candidates is divided into two parts. First, a loose selection is performed 
that retains the majority of signal events whilst reducing the background substantially. Af¬ 
ter this the —>■ 'ip{2S)K^TT~ peak is clearly visible. Subsequently, a multivariate method 
is used to further improve the signal-to-background ratio and to allow the observation of 
the 5° —)■ '0(25')iF+7r“ decay. 

The selection starts by reconstructing the dimuon decay of the 0(25') meson. Pairs 
of oppositely charged particles identihed as muons with px > 550MeV/c are combined to 
form 0(25') candidates. The invariant mass of the dimuon pair is required to be within 
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To form candidates, the selected '0(25') 


60MeV/c^ of the known 0(25) mass 
mesons are combined with oppositely charged kaon and pion candidates. Tracks that do 
not correspond to actual trajectories of charged particles are suppressed by requiring that 
they have px > 250MeV/c and by selecting on the output of a neural network trained 
to discriminate between these and genuine tracks associated to particles. Combinatorial 
background from hadrons originating in the primary vertex (PV) is suppressed by requiring 
that both hadrons are signihcantly displaced from any PV. Well-identihed hadrons are 
selected using the information provided by the Cherenkov detectors. This is combined 
with kinematic information using a neural network to provide a probability that a particle 
is a kaon (V^), pion {V") or proton {V^). It is required that is larger than 0.1 for the 
candidate and that is larger than 0.2 for the vr” candidate. 

A kinematical vertex fit is applied to the candidates ^ . To improve the invariant 
mass resolution, the fit is performed with the requirement that the candidate points 
to the PV and the ‘ip{2S) is mass constrained to the known value 26 . A good quality of 


the vertex ht Xdtf: is required. To ensure good separation between the and B^ 
signals, the uncertainty on the reconstructed mass returned by the ht must be less than 
11 MeV/c^. Combinatorial background from particles produced in the primary vertex is 
further reduced by requiring the decay time of the meson to exceed 0.3 ps. 

Four criteria are applied to reduce background from specihc 6-hadron decay modes. 
First, the candidate is rejected if the invariant mass of the hadron pair calculated assuming 
that both particles are kaons is within lOMeV/c^ of the known 0 meson mass 26 , 
suppressing B^ —>■ 0(25)0 decays where one of the kaons is misidentihed as a pion. Second, 
to suppress B^ —>■ 0(25)7r’''7r“ events where one of the pions is incorrectly identihed as 
a kaon, it is required that for the kaon candidate. This rejects 80 % of the 

background from this source whilst retaining 90% of signal candidates. Third, to 
suppress background from —)■ 0(25)p7r“ decays where the proton is misidentifed as 
a kaon, candidates with > 0.3 and an invariant mass within 15MeV/c^ of the known 
A? mass 


26 


are discarded. Finally, to reduce background from a B^ —)■ ip(2S)K^ decay 
combined with a random pion, candidates where the reconstructed 'ip{2S)K^ invariant 
mass is within 16MeV/c^ of the known B^ mass 26 are rejected. Background from the 


decay A° —)■ 0(25)piF“ with misidentihed hadrons does not peak at the 5° mass and is 
modelled in the ht. 

To further improve the signal-to-background ratio, a multivariate analysis based on 
a neural network is used. This is trained using simulated B^ signal events together 
with candidates from data with a mass between 5500 and 5600MeV/c^ that are not used 
for subsequent analysis. Eight variables that give good separation between signal and 
background are used: the number of clusters in the large-area silicon tracker upstream 
of the magnet, for the kaon candidate, V" for the pion candidate, the transverse 
momentum of the B^^^y the minimum impact parameter to any primary vertex for each of 
the two hadrons, Xdtf the hight distance in the laboratory frame of the B^y candidate 
divided by its uncertainty. The ratio Ng/y/Ns + is used as a hgure of merit, where 
Ns (Nb) is the number of signal (background) events determined from the invariant mass 
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fit (see Sect. |^. The maximum value of this ratio is found for a threshold on the neural 
network output that rejects 98% of the background and retains 81% of the signal for 
subsequent analysis. 


4 Invariant mass fit 


A maximum likelihood ht is made to the unbinned ijj{2S)K~^7r invariant mass distribution, 
m{'ijj{2S)K~^7r~), to extract the and signal yields. The B^ signal component is 


modelled by the sum of two Crystal Ball functions 28 with common tail parameters and 


an additional Gaussian component, all with a common mean. All parameters are hxed 
to values determined from the simulation apart from the common mean and an overall 
resolution scale factor. The simulation is tuned to match the invariant mass resolution 
seen in data for the B~^ —>■ J/iIjK^ and B^ —)■ J/'il)K~^7r~ decay modes. Consequently, the 
resolution scale factor is consistent with unity in the ht to data. The B^ component is 
modelled with the same function, with the mean value of the B^ meson mass left free in the 
ht. The resolution parameters in this case are multiplied by a factor of 1.06, determined 
from simulation, which accounts for the additional energy release in this decay. 

The dominant background is combinatorial and modelled by an exponential function. 
A signihcant component from B^ —)■ 'ip{2S)(j) decays is visible at lower masses than 
the B^ peak. This is modelled in the ht by a bifurcated Gaussian function where the 
shape parameters are constrained to the values obtained in the simulation and the yield 
constrained to the value determined in data under the hypothesis that both hadrons are 
kaons. Additional small components from —>■ 'iP{2S)tt^tt~ and —>■ 'ip{2S)pK~ decays 

are modelled by bifurcated Gaussian functions. The shapes of these components are hxed 
using the simulation and the yields are determined by normalising the simulation samples 
to the number of candidates for each modes found in data using dedicated selections. 
Contributions from partially reconstructed decays are accounted for in the combinatorial 
background. In total, the ht has ten free parameters. Variations of this ht model are 
considered as systematic uncertainties. 

Figure]^ shows the invariant mass distribution observed in the data together with the 
result of a ht to the model described above. Binning the data, a y^-probability of 0.30 
is found. The moderate mismodelling of the B^ peak is accounted for in the systematic 
uncertainties. The ht determines that there are 329 ± 22 decays and 24207 ± 160 B^ 
decays. The B^ —>■ 'il>{2S)K~^7r~ mode is observed with high signihcance. 

The precision of the momentum scale calibration of 0.03% translates to an uncertainty 
on the B^ and B^ meson masses of 0.3MeV/c^. Therefore, it is chosen to quote only the 
mass diherence in which this uncertainty largely cancels, 

M(B°) - M(5°) = 87.45 ± 0.44 (stat) ± 0.09 (syst) MeV/c^ 

This procedure has been checked using the simulation, which gives the input mass difference 
with a bias of 0.05MeV/c^ that is assigned as a systematic uncertainty. Further systematic 
uncertainties arise from the momentum scale and mass ht model. Varying the momentum 
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Figure 2: Invariant mass distribution for selected 'ip{2S)K'^7r~ candidates in the data. A fit to 
the model described in the text is superimposed. The full fit model is shown by the solid (red) 
line, the combinatorial background by the solid (yellow) and the sum of background from the 
exclusive b —)• iIj{2S)X modes considered in the text by the shaded (blue) area. The maximum of 
the y-scale is restricted so as to be able to see more clearly the —)• ^j;{2S)K~^^T~ signal. The 
lower plot shows the differences between the fit and measured values divided by the corresponding 
uncertainty of the measured value, the so-called pull distribution. 


scale by 0.03% leads to an uncertainty of 0.04MeV/c^. The effect of the fit model is 
evaluated by considering several variations: the relative fraction of the two Crystal Ball 
functions is left free; the slope of the combinatorial background is constrained using 
candidates where the kaon and pion have the same charge; the Gaussian constraints 
on the background from the —)■ 'ip{2S)(j) mode are removed; and the tail parameters 
of the Crystal Ball functions are left free. The largest variation in the mass splitting 
is 0.06 MeV/c^. The total systematic uncertainty is given by summing the individual 
components in quadrature. 


5 Amplitude analysis 

Figure shows the Dalitz plot of the selected —)■ ip(2S)K^'K~ candidates in the signal 
range, m{'ip{2S)K^'K~) G [5350,5380] MeV/c^. There is a clear enhancement around the 
known K*{892)^ mass and no other significant enhancements elsewhere. To determine 
the fraction of decays that proceed via the iC*(892)° resonance, an amplitude analysis 
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Figure 3: Dalitz plot for the selected —)• 'ip{2S)K~^Tr candidates in the signal window 

m{'il^{2S)K+7r-) E [5350, 5380] MeV/c^. 



Figure 4: Definition of the helicity angles. 


is performed, similar to that used in Ref. for the analysis of the —>■ ilj{2S)K~^7r~ 
mode. The hnal-state particles are described using three angles Q = (cos 6*^^, cos 6*^, 0) 
in the helicity basis, defined in Fig. and the invariant K~^n~ mass, mx-K = m{K^Ti~). 
The total amplitude is S{mK-K-i^)s{jnK'K-i^) + B{mKn,^), where S{mKn,^) represents 
the coherent sum over the helicity amplitudes for each considered resonance or 

non-resonant component. The detection efficiency, e{mKTT,cos6x,cos6^,(j)), is evalnated 
using simulation and parameterised using a combination of Legendre polynomials and 
spherical harmonic moments, given by 

£{mK-K, cos Ok, cos e ^, 41 ) = ^ c'^’''"^Pa{cos9K)Ybc{cos9^,(j))Pd ~ > (1) 

a,b,c,d ^ ^Kn / 

where is the minimum (maximum) value allowed for rriKiT in the available phase 

space of the decay. The coefficients of the efficiency parameterisation are computed by 
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Figure 5: Distributions of (a) cos0^, (b) (/)•, (c) cosOk and (d) m{K~^7r~) of simulated —)• 
^lJ{2S)K~^^T~ decays in a phase space configuration (black points) with the parameterisation of 
the efficiency overlaid (blue lines). 


summing over the Nmc events simulated uniformly in the phase-space as 


^ahcd 


Nj 


MC 


^ 2a + 1 2d + 1 ^ 
9 9 


COS 9Ki)Ybc{cos e^., (j)i)Pd 




m 


Ktt , 


Ktt 


m 


Ktt 



C 

9i' 

( 2 ) 


where Qi = piQi, with pi {qi) being the momentum of the K^'k~ system meson) in 
the (ib+vr”) rest frame and C* is a normalising constant with units GeV^/c^. This 
approach provides a description of the multidimensional correlations without assuming 
factorisation. In practice, the sum is over a finite number of moments (a < 2, 6 < 2, c < 2 
and d <2) and only coefficients with a statistical signihcance larger than hve standard 
deviations from zero are retained. The one-dimensional projections of the parameterised 
efficiency are shown in Fig. superimposed on the simulated event distributions. 

The background probability density function, is determined using a similar 

method as for the efficiency parameterisation. In this case the sum in Eq. ([^ is over the 
selected events with m{'ip{2S)K~^n~) > 5390MeV/c^ and Qi = 1. Only moments with a < 2, 
6 = 0, c = 0 and d <2 and a statistical signihcance larger than hve standard deviations 
from zero are retained. The one-dimensional projections of the parameterised background 
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Figure 6: Distributions of (a) cos0^, (b) (j), (c) cos9k and (d) m{K~^7r ) of —>• 'ijj{2S)K~^7T 

candidates with m{'ijj{2S)K~^7r~) > 5390MeV/c^ (black points), with the parameterisation of the 
background distribution overlaid (blue lines). 


distribution are shown in Fig. superimposed on the sideband data. As a consistency 
check, the {mKn, ^) distributions for events with m{'ip{2S)K~^7i~) > 5390MeV/c^ are found 
to be compatible with the same distributions obtained from a like-sign {'il>{2S)K^n^) 
sample. 

The default amplitude model is constructed using contributions from the A*(892)'^ 
resonance and a K^tt~ S-wave modelled using the LASS parameterisation j^. The 
magnitudes and phases of all components are measured relative to those of the zero helicity 
state of the A*(892)° meson and the masses and widths of the resonances are hxed to 
their known values |^. The remaining eight free parameters are determined using a 
maximum likelihood fit of the amplitude to the data in the signal window. The background 
fraction is fixed to 0.28, as determined from the fit described in Sect. The fit fraction for 
any resonance R is defined in the full phase space, as fn = J Sr drriKndSl / J S dmKndSl, 
where Sr is the signal amplitude with all amplitude terms set to zero except those for R. 
The fractions of each component determined by the fit are /ir*( 892 )o = 0.645 ± 0.049, and 
/s-wave = 0.339 ±0.052, where the uncertainty is statistical only. The fractions do not sum 
to unity due to interference between the different components. Variations of the S-wave 
description and default mixture of A+vr” resonances, including the introduction of the 
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Figure 7: Distributions of (a) cos9^, (b) (j), (c) cosOk and (d) m{K^'K~) for selected —)• 
'iI)[2S)K'^'k~ candidates (black points) with the projections of the fitted amplitude model overlaid. 
The following components are included in the model: iF*(892)^ (red dashed), LASS S-wave 
(green dotted), and background (grey dashed-dotted). The residual pulls are shown below each 
distribution. 


spin-2 7 ^ 2 ( 1430 )° meson or an exotic Z~ meson, are considered but found to give larger 
values of the Poisson likelihood degree of freedom or lead to components with ht 

fractions that are consistent with zero. For each model the number of degrees of freedom is 
calibrated using simulated experiments. The variations in amplitude model are considered 
as sources of systematic uncertainty. The longitudinal polarisation fraction of the iF*(892)° 
meson is dehned as /l = HU {Hi + Hl_ + Hi), where are the magnitudes of the 

K*(892)^ helicity amplitudes. This is measured to be /l = 0.524 ± 0.056, where the 
uncertainty is statistical. The projections of the default ht for the helicity angles and 
invariant K~^n~ mass are shown in Fig. 

5.1 Systematic uncertainties of amplitude analysis 

A summary of possible sources of systematic uncertainties that affect the amplitude 
analysis is reported in Table The size of each contribution is determined using a set of 
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Table 1: Summary of systematic uncertainties on the itr*(892)*^ fit fraction and /l- Rows marked 
with (*) refer to uncorrelated sources of uncertainty between the and modes for the 
computation of the ratio of branching fractions. 


Source K*{892)^ fit fraction /l 


(*) A+vr” amplitnde model 

0.028 

0.017 

(*) S-wave model 

0.018 

0.010 

K* resonance widths 

0.005 

0.008 

Blatt-Weisskopf radius 

0.014 

0.003 

Breit-Wigner parameters (m^ vs. rriK-K) 

0.026 

0.005 

(*) Background parameterisation 

0.014 

0.012 

(*) Background normalisation 

0.007 

0.011 

Efficiency model (parameterisation) 

0.011 

0.007 

Efficiency model (nenral net) 

0.002 

0.004 

Quadrature sum of systematic uncertainties 

0.049 

0.029 

Quadrature sum of uncorrelated systematic nncertainties 

0.037 

0.026 

Statistical nncertainty 

0.049 

0.056 


simulated experiments, of the same size as the data, generated under the hypothesis of 
an alternative amplitude model. These are htted once with the default model and again 
with the alternative model. The experiment-by-experiment difference in the measured £t 
fractions and /l is then computed and the sum in qnadrature of the mean and standard 
deviation is assigned as a systematic nncertainty to the corresponding parameter. 

The systematic dependence on the K^'k~ amplitnde model is determined using the 
above procedure, where the alternative model also contains a spin-2 iT2(1430) component. 
This leads to the dominant systematic uncertainty on the it'*(892)° ht fraction and /l. 
The systematic dependence on the K^tt~ S-wave model is determined nsing simnlated 
experiments where a combination of a non-resonant term and a contribution is 

used in place of the LASS parameterisation. In addition, the amplitnde model contains 
parameters that are hxed in the default fit such as the masses and widths of the resonances 
and the Blatt-Weisskopf radins. The radins controls the effective hadron size and is set 
to 1.6 (GeV/c)“^ by defanlt. Alternative models are considered where this is changed to 
3.0 (GeV/c)-i and 0.8 (GeV/c)-L 

A large source of systematic uncertainty comes from the choice of convention for the 
mass, m, in the {p/m)^^ terms of the amplitude. The default amplitude model follows the 
convention in Ref. by using the resonance mass. This is different to that in Ref. |[^ 
where the rnnning resonance mass (rriKn) is used in the denominator. This choice is 
motivated by the improved fit qnality obtained when using the resonance mass. 

The systematic uncertainty related to the combinatorial background parameterisation 
is determined using an amplitude model with an alternative background description that 
allows for higher moment contributions (a < 2, 6 < 2, c < 2 and d < 2). The combinatorial 
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background normalisation is determined from the fit to the m{'ip{2S)K^T[~) distribution 
and is hxed in the amplitude ht. The systematic uncertainty related to the level of 
the background is estimated by using an amplitude model with the background fraction 
modihed by ±10%. 

The efhciency parameterisation is tested by re-evaluating the coefficients, allowing for 
higher order moments (a < 4, 6 < 4, c < 4 and d < 4). Similarly, to test the dependence of 
the efficiency model on the neural network requirement, an alternative model is used with 
the efficiency parameterisation determined from the simulated events that are selected 
without applying the requirement. There is a negligible systematic uncertainty caused by 
the lifetime difference between the and mesons. 


6 Branching fraction results 


Two ratios of branching fractions are calculated, B{bI —>• ^jJ{2S)K^^l )/B{B^ —)■ 
i:{2S)K+7r-) and B(bI {2S)K*(892)^)/B{B^ i;{2S)K*(892)^). These are deter¬ 

mined from the signal yields given in Sect, [^correcting for the relative detector acceptance 
using simulation. The simulated B^ samples are reweighted with the results of the angular 
analysis presented in Sect. [^ Similarly, the B^ simulated data are reweighted to match 
the results given in Ref. j^. For the inclusive branching ratio, the relative efficiency 
between the two modes is found to be 0.975 ±0.014 whilst for the K*{892)^ component it is 
1.027 ± 0.021. The uncertainty on these values is propagated to the systematic uncertainty. 

Since the same hnal state is considered in the signal and normalisation mode, most 
sources of systematic uncertainty cancel in the ratio. The remaining sources are discussed 
in the following. The variations of the invariant mass ht model described in Sect, j^are 
considered. The largest change in the ratio of yields observed in these tests is 3.7%, which 
is assigned as a systematic uncertainty. Differences in the px spectra of the B^ meson are 
seen comparing data and the reweighted simulation. If the px spectrum in the simulation 
is further reweighted to match the data, the efficiency ratio changes by 0.7%, which is 
assigned as a systematic uncertainty. 

To test the impact of the chosen K^tt~ amplitude model for the B^ channel, the 
simulated events are reweighted using a model consisting of the K*{892)^ resonance, 
the LASS description of the S-wave and the 772(1430) resonance. This changes the 
efficiency ratio by 0.6%, which is assigned as a systematic uncertainty. To calculate the 
77* (892)° branching ratio, the fraction of candidates from this source is needed. For the 7?° 
channel this is determined from the amplitude analysis to be 0.645 ± 0.049 ± 0.049 and the 
corresponding fraction for the 7?° channel is 0.591 ± 0.009 [^, leading to a 6.0% systematic 
uncertainty. All of the uncertainties discussed above are summarised in Table The 
limited knowledge of the fragmentation fractions, fs/fd = 0.259 ± 0.015 
an uncertainty of 5.8%, which is quoted separately from the others. 


31-33 , results in 
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Table 2: Systematic uncertainties on the ratio of branching fractions. 


Relative uncertainty % 


Source 

Inclusive 

77* (892)° 

Simulation sample size 

1.4 

2.2 

Fit model 

3.7 

3.7 

Detector acceptance 

0.7 

0.7 

amplitude model 

0.6 

- 

77* (892)° fit fraction 

- 

6.0 

Quadrature sum 

4.1 

7.4 


7 Summary 


Using a data set corresponding to an integrated Inminosity of 3.0 fb ^ collected in pp 
collisions at centre-of-mass energies of 7 and 8 TeV, the decay —)■ 'ip{2S)K^7i~ is 

observed. The mass splitting between the and B^ mesons is measured to be 

M(B°) - M(5°) = 87.45 ± 0.44 (stat) ± 0.09 (syst) MeV/cl 


This is consistent with, though less precise than, the value 87.21 ± 0.31MeV/c^ obtained 
by averaging the results in Refs. [^|^. Averaging the two numbers gives 

- M(R°) = 87.29 ± 0.26MeV/cl 


The ratio of branching fractions between the B^ and B^ modes is measured to be 


B{bI ^lJ{2S)K+7l-) 
13{B0 ij{2S)K+7r-) 


5.38 ± 0.36 (stat) ± 0.22 (syst) ± 0.31 {fs/fd) %• 


The fraction of decays proceeding via an intermediate 77*(892)° meson is measured with 
an amplitude analysis to be 0.645 ± 0.049 (stat) ± 0.049 (syst). No significant structure is 
seen in the distribution of m{'ip{2S)7i~). 

The longitudinal polarisation fraction, /l, of the 77*(892)° meson is determined 
as 0.524 ± 0.056 (stat) ± 0.029 (syst). This is consistent with the value measured 
in the corresponding decay that proceeds through an intermediate J/'ijj meson, 
/l = 0.50 ± 0.08 (stat) ± 0.02 (syst) j^. The present data set does not allow a test of 
the prediction given in Ref. that /l should be lower for decays closer to the kinematic 
endpoint. 

Using the 77*(892)° fraction determined in this analysis for the 7?° component, the 
corresponding number for the 7?° mode from Ref. j^|, and the efficiency ratio given in 
Sect. the following ratio of branching fractions is measured 


B{bI ^/>(2R)77*(892)°) 
i3(R° ^^/>(2R)i7*(892)°) 


5.58 ± 0.57 (stat) ± 0.40 (syst) ± 0.32 {fs/fd) %• 
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The -Bg —)■ iI){2S)K^t:~ mode may be useful for future studies that attempt to control the 
size of loop-mediated processes that influence CP violation studies and offers promising 
opportnnities in the search for exotic resonances. 
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